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We investigate the statistics of the gap, G„, between the two rightmost positions of a Markovian 
one-dimensional random walker (RW) after n time steps and of the duration, L n , which separates 
the occurrence of these two extremal positions. The distribution of the jumps ryi's of the RW, f{rj), 
is symmetric and its Fourier transform has the small k behavior 1 — f(k) ~ |fc| M with < n < 2. 
We compute the joint probability density function (pdf) P n (g,l) of G n and L n and show that, 
when n — > oo, it approaches a limiting pdf p(g,l). The corresponding marginal pdf of the gap, 
Pgap(s)) is found to behave like p gap (g) ~ 3 _1_M for j > 1 and < /i < 2. We show that the 
limiting marginal distribution of L„, ptime(O) has an algebraic tail ptime(0 ~ ?~ 7 ' M) for / 3> 1 with 
7(1 < n < 2) = 1 + and 7(0 < jj, < 1) = 2. For I, g » 1 with fixed lg~^, p{g, I) takes the scaling 



form p(g, I) 



p^ilg M ) where is a (/^-dependent) scaling function. We also present 



numerical simulations which verify our analytic results. 
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Extreme and order statistics are currently the sub- 
ject of numerous studies in various areas of sciences. In 
many circumstances, statistical systems are not governed 
by typical or average events but instead by anomalously 
rare and intense ones Illustrative examples of such 
situations include for instance natural disasters or fi- 
nancial crisis Q where extreme events, like earthquakes, 
tsunamis or financial crashes may have drastic conse- 
quences. It is also now well established that extreme 
value questions play an important role in the statistical 
mechanics of disordered systems ■ 

When studying extreme statistics, one is usually inter- 
ested in studying the maximum A max among a collection 
of N + 1 random variables Xo,X\, ■ ■ ■ ,Xn. However, 
the knowledge of the statistics of this single global vari- 
able, though important, does not always provide enough 
valuable information. In particular, a crucial question 
concerns the crowding of events near the maximum. For 
instance, a rare event like an earthquake is usually not 
isolated but is followed (or preceded) by smaller ones, 
called aftershocks (respectively foreshocks) [H, Q . Fore- 
shocks and aftershocks are also known to occur before 
and after a financial market shock [l^, [ll|. This is a 
natural question in statistical physics too, when one is 
interested not only in the ground state properties of a 
disordered system but also in the finite low temperature 
physics, which involves the low lying energy states, close 
to the ground state [IH- This has led to the study of 
the density of near extreme events, both in statistics [l3j 
and in physics [l4| . which essentially counts the number 
of events Xi 's which are at given distance from X max 14 1 . 



nomenon is to study the order statistics of the Xi's, that 
is the statistics of X max = Mi >n > ... > Mk, n > ■■■ > 
Mn+x^n where Mk, n denotes the k th maximum of the 
set {Xq, Xi, . . . , X n }. One natural question for disor- 
dered systems, if one interprets — Mfc : „ as the k th energy 
level of the system, is the distribution of the first gap 
G n = M\ t n — M2, n as it controls, to a large extent, the 
low temperature properties of the system. The gap G n is 
also an important quantity for applications in seismology 
as it can model the difference in magnitude between the 
mainshock and its largest aftershock or foreshock [H-Q1]. 
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Crowding of near extreme events has also been studied in 
the context of sporting events, like marathon packs [lij . 

Another natural way to characterize this crowding phc- 



FIG. 1. Realization of a random walk JT]) of n steps. Here we 
focus on the joint pdf P n (g,l) of G n and L n in the limit of 
large n. 



Aside from the gap, another key random variable is the 
time L n elapsed between these two extreme events (see 
Fig. []}, which is particularly important for the statistics 
of earthquakes or financial crashes. 

While the study of the fluctuations of G n and L n is 
very well understood in the case of independent and 
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identically distributed (i.i.d.) random variables [20[, this 
question is highly non trivial for strongly correlated vari- 
ables Xi's. Yet it is known, for instance, that after- 
shocks exhibit long range correlations [1, B (both 
spatial and temporal) , in which case a model of indepen- 
dent variables cannot predict anything sensible about the 
distribution of the time separating the mainshock from 
the largest aftershock. The importance of order statis- 
tics for strongly correlated variables came up recently in 
several other physical contexts, notably in the study of 
the branchin g Br ownian motion (BBM) [23| and also for 
signals [24j| . In the former case, G n represents the 
gap between the two rightmost particles. The pdf of G n 
was studied in [23I ] , though an exact analytical expression 
remains a hard task. Therefore any exact results for the 
statistics of the gap and the time between the two first 
maxima for a set of strongly correlated variables would 
be highly desirable. 

In this Letter, we obtain exact results and find a very 
rich behavior for the statistics of G n and L n in the case 
where the X.^s correspond to the positions of a ran- 
dom walker (RW) at discrete times i's. Such RW is cer- 
tainly the simplest, yet non trivial, set of strongly cor- 
related random variables for which interesting problems 
of records j2f| and order statistics [26| can be solved ex- 
actly. Hence, while RW might not be a realistic model 
for earthquakes, this is a useful laboratory where the ef- 
fects of correlations on the statistics of G n and L n can 
be studied in detail. 

In our model, the RW starts at Xq = at time and 
evolves via 



Xi — X; 



Vi 



(1) 



where the 77.^ 's are i.i.d. random jumps each drawn from 
a symmetric distribution f(rj). Its Fourier transform, 
f(k) = e %kn f(r])dr], has the small k behavior 



f(k) = 1 - |ofc|* + odfcl") , < n < 2 , 



(2) 



where a is the characteristic length scale of the jumps. 
In particular, for < /1 < 2, one has f(rj) ~ B/lrj]^ 1 , 
for large \r)\, with B = sin (fiir/2)T\pL + 1]/tt. Let Mi,„ 
and M2 >n respectively denote the first and second max- 
ima of the random walk ([T]) after n steps and write n\ 
and ri2 the times at which they are reached: X ni = M\^ n 
and X n2 = M2, n - The purpose of this work is to study 
the joint pdf P n (g, /) of the gap G n = M\ >n — M 2j „ and 
the time L n = n\ — ni between the occurrence of these 
first two maxima (Fig. [TJ. 

It is useful to summarize our main results. We first 
show that P n (g, has a well defined limiting pdf p(g, I) 
as n — > 00. By integration over / one obtains an exact 
expression for the marginal distribution p gap (g) <[T7| SH I 
the full form of which depends on f(rj). For < \i < 2, 
we show that p gap (g) has an algebraic tail 



Pgap(ff) ~ Bp g 



-l-n 



g -> 00 , 



(3) 



where B^ is a computable constant. For fi = 2, the full 
distribution can be computed in some specific cases only, 
and the tail itself remains non-universal and sensitive to 
f(rj). By integration over g, one finds that the marginal 
distribution ptime(0 displays an algebraic tail whose ex- 
ponent depends only on the Levy index [i as follows: 



T^timc 



Ail- 1 - 1 ^ , 1< fi < 2 
Au (log/)/- 2 , u = l 



(0~ <An(\ogl)l- 2 



, A m I 



, 0< M ^i<l, 



(4) 



for / — > 00, where the amplitudes Ai, An and Am are non 
universal and k £ N [23] ■ The third line of (@| reveals an 
unexpected freezing phenomenon of the exponent char- 
acterizing the algebraic tail of ptimc(Z) as /i decreases past 
the value /i c = 1. Interestingly, we see that the first mo- 
ment of j?time(Z) is not defined. This means that although 
the typical size of L n is 0(1), its average diverges with 
n. From (j4|) one can estimate that (L n ) ~ n 1-1 ^ for 
1 < fi < 2, while (L n ) ~ logn for < fi < 1. Finally, 
in the scaling regime j, / > 1 with fixed lg~^ and for 
< /X < 2, we find the following scaling form 



p(g,i) ~ g-^P^lg-*) 



(5) 



where the //-dependent scaling function Pn{y) is inte- 
grable over [0, +00) with the algebraic tail 



Py.{y) ~ V 



-1-1/ n 



y 



(6) 



The starting point of our analysis is an exact formula 
for the joint pdf of the random variables ni, i%2, and 
G n = M\, n — M2 t n f° r a random walk of n steps, in terms 
of the following two central objects. The first one is the 
survival probability q n (x) for a random walker, starting 
at x > 0, to stay on the positive axis up to step n. Note 
that by q n {x) = Pr(Mi „ < x), it is also the cumulative 
distribution of Mi_ n [28j |, A complete characterization of 
q n {x) is given by the Laplace transform (LT) with resp ect 
to (wrt) x of its generating function (GF) wrt n 
(see also (HI), 



respect 

la El 



/ dx V q m (x)s m e- X * = <p(s, A) 

Jo „^ XVT^l 



where the function y>(s, A) is given by 



tp(s,X) 



exp -- 

\ 7T 



A r°°ln[l-sf(k)] 



k 2 + X 2 



dk 



(7) 



(8) 



The second object is the probability p n (x)dx for a ran- 
dom walker, starting at x — 0, and conditioned to stay 
positive, to arrive at step n in the interval [x,x + dx]. 
The counterpart of for p„ (x) reads [Hf (see also (33l |) 



POO 

/ dxY,PMs m e~ Xx 



tp(s,X) . 



(9) 
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To compute P(g, n%, n 2 \n), we exploit the renewal prop- 
erty of the random walk and divide the interval [0, n] into 
three independent parts (see Fig. [T]) of duration li = n 2 , 
h = n i ~ n 2 and I3 = n — n 1; (here we suppose without 
loss of generality that n\ > n 2 ). One has, 

P(g, n u n 2 \n) = w 1 (li)w 2 {l)w 3 (l 3 )Si 1+ i + i atn , (10) 

where Wk{lk), with k = 1,2,3 and l 2 = I, denotes the 
weight of the paths in each of the three subintervals (see 
Fig. [1]) and 5ij is the Kronecker delta function. The 
weight wi(li) is simply given by the survival probability 

«;i(/i) = g Ii (0). (11) 

This can be readily seen by reversing the direction of 
both space and time axis and taking as a new origin the 
point of coordinates (h,M 2 ^ n ) (notice that P(g, fi\, n 2 \n) 
is obtained by integrating over the value of M 2tU keeping 
the value of the gap g fixed). To compute w 2 (l) we isolate 
the last step, of amplitude g + y, before the maximum 
Mi^ n is reached, from the first n 2 — n\ — 1 steps on this 
interval (Fig. [IJ. The weight associated to these n 2 — 
ni — 1 steps is given by p n2 -n 1 -i{y) and hence 

POO 

u*(i)= / P i-i(y)f(g + y)dy . (12) 
Jo 

Similarly, to compute the weight ^3(^3) we isolate the 
first step, of amplitude g + x, after the maximum M\ >n is 
reached from the last n — n\ — \ steps. The weight asso- 
ciated to these last steps is simply given by q n — ni -\(x) 
and one has 

/>oo 

w 3 (h) = q l Ax)f{g + x)dx . (13) 
Jo 

The joint pdf P n (g,l) is obtained formally from (fT0|) 
as P n (g,l) = T,t7s=o P (9^ n ^ n 2\n), and using (JTTJ) - (JTHJ 
together with (|7|) and ©, one obtains an explicit expres- 
sion for the double GF of P n (g,l) wrt n and I [we recall 
that P n (g,l) = P n (g, — I)}. Namely, 

oc n 2 poo 

s n J2 fl P ^9, l) = — / u(st, y)f(g + y) dy 

poo 

x / h(s,x)f(g + x)dx , (14) 
Jo 

where u(t,y) and h(s,x) are the inverse LTs of (p(t,X) 
and cp(s, A)/A, respectively: 

f™u(t,y)e- x ydy = ip(t,\) , 

(15) 

J °° h(s,x)e Xx dx = ip(s, A)/A . 

From the s — > 1 limit of (fT4|) one can extract the large 
n limit of P n (g,l). It is readily seen that the leading 
divergence on the right hand side of (fl~4"l) is a simple pole 




FIG. 2. a: Plot of p g a P ((?) for f(rf) = ^\x\ exp ( — \x\) and 
n = 10 4 . The dotted line corresponds to the exact result 
given in the text given below (|18l) . b): Plot of Pga P (<?) for 
pi = 1.4 and n = 10 3 . The solid line is the exact asymptotic 
result. 

at s = 1: this implies that P n (g, I) converges to a limiting 
distribution p(g, I) as n — > 00 with GF wrt / 

00 

p(g, z) = J2 zl P(9, - h(z, g)h(g) (16) 
i=o 

where 

h(z,g) = z u(z,y)f(g + y)dy , 

(17) 

h{g) = / °° h(l, x)f(g + x) dx . 

Expression (fl~6l) . together with (fl~7|) . is the central result 
of our study from which the various behaviors announced 
in the introduction can be derived. 

We first focus on the marginal distribution of the gap 
Pgap(<?) which, for any jump distribution, is exactly given 
by 

Psap (g)=2p(g,l) = 2I 1 (l,g)I 2 (g), (18) 

where the factor of 2 comes from the configurations with 
/ 0. For /i = 2, there are some particular cases in which 
(fT5|) can be computed explicitly. For instance, if f(rj) = 
(b/2) exp (— &I77I) one finds p gap {g) = 26exp(— 2bg), and 
for/(r?) = (6 2 /2)|77|exp(-6|77|)onehasp gap (g) = c[(V3+ 
2bg) 2 - l}e~ 2b 3 with c = 2b/ (1 + VS) 2 . Fig. la) shows 
a numerical check of the latter exact result. On these 
two examples we see that, for pL = 2, even the tail of 
Pg a p(g) depends on the the details of the jump distri- 
bution. On the other hand, for < \i < 2, the tail of 
Pg a p(g) depends on the Levy index pi only. In the large 
g limit, it turns out that the integrals over x and y in 
Eq. ([!"?)) are dominated by large values of x,y ~ O(g). 
Hence, to study the large g limit, one needs the large 
argument behavior of u{l,y) and h{\,x) in Eq. (|15[) . 
These behaviors can in turn be obtained by analyzing 
the small A behavior of <p(l,A) in Eq. (jSJ. One finds 
that u{l,y) - bxy^' 2 - 1 , with h = [a^ 2 T(pi/2)}~\ and 
h(l,x) - ax^' 2 with ci = K/ 2 r(l + m/2)] -1 - Using 
these asymptotic behaviors, it can be shown from (|17|) 
and (jTSJ) that p E&v {g) ~ B^g" 1 "^ as announced in Eq. 
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© with the amplitude = a^>/[r(l - /i/2)] 2 . Fig. 
([2]) b) shows a numerical estimate of p gap (g) for /i = 1.4: 
the tail behavior is in good agreement with the analytical 
predictions 

We now come to the marginal distribution ptime(0 of 
the time I elapsed between the first two maxima. Its GF 
is readily obtained by integrating (|16[) over g. One finds, 

Pth-nc(z) = ^2 ^Pthne (0 = / dgh(z, g)h(g) ■ (19) 

(=0 J ° 

For an exponential jump distribution, f(rf) = 
(b/2) exp (— b\r]\), corresponding to fi = 2, Ptime(0 can 
be computed exactly, Ptime(0 = I\Z — ^)/[4y / 7tT(/ + 1)] 
and ptimo(0 oc l~ 3 ^ 2 , for large I. For other jump distri- 
butions, an explicit computation of Ptime(0 is generally 
impossible but its large I behavior can be obtained from 
the z — > 1 limit of the GF (fT9"|). This analysis is rather 
subtle and requires the analysis of tp(z,\) when z — > 1. 
In this limit, Eq. ((8} yields 

rtz,\) = <p(l,\) + a„^^ l p(l,\) + o^^} , 

(20) 

with a M = (a sin (tv//j)) . For 1 < /i < 2, the asymp- 
totics can be directly used to compute the large I 
behavior of ptimc(0- One finds the first line of ((4]) with 

poo 

Ai = (a7r) _1 r(l + l//x) / [7 2 (.g)] 2 dff. (21) 
Jo 

A numerical check of these results for 1 < fj, < 2 is shown 
in Fig. H For /i < 1 (and fi^l/k with fc e N and fc > 1), 
things are more complicated since the (/-integral defining 
Ai diverges (because ^2(5) ~ .9 for large g). In this 




FIG. 3. Left: Log-log plot of f>time(0 for different values of fi 
and n = 10 (the case [i = 2 corresponding to an exponential 
jump distribution). The straight lines correspond to an alge- 
braic decay as predicted in Eq. Q (note that 1 + 1/1.4 ~ 1.7). 
These data have been vertically shifted for clarity. 

regime, the large / behavior of ptime(0 can be obtained 



by the following scaling argument. The expansion (f20|) 
holds for (1 — z)~ j\ <C 1. This suggests that the g- 
integral in (f2~Tj) should actually been cut-off around l 1 ^ 
such that for [1 < 1, the second line of (U]) is replaced 
with ptimc(0 ~ Aml~ 2 , hence the freezing of the tail of 
Ptime(0 for H < 1, as discussed below Eq. (|4]). For the 
marginal case fj, = 1, one finds the logarithmic correction 
given in the second line of (j4]). Fig. [3] shows numerical 
results for ptime(0 with /i = 0.4 and /1 = 0.3 which both 
corroborate ptime(0 oc l~ 2 . 

From Eqs. ([3]) and Q it is possible to get the scal- 
ing form of the joint pdf p(g, I) for large g and / in ([5]). 
According to standard scaling arguments, p(g, I) is ex- 
pected to depend on / through the dimcnsionlcss com- 
bination lg~^. Moreover, for large g we showed that 
Pgap(g) ~ g 1 ^- From these two arguments, it is natu- 
ral to expect that p(g, I) has the scaling form ([5]) where 
the function p^{y) is integrable over [0, +00). Indeed, 
one can easily check that integrating ([5]) over I yields 
P gap (5) ~ .g" 1 "^ / °° Pn(y)dy, in agreement with the large 
g behavior p gap (g) ~ g~ 1 ~ t * in Eq. , as it should. 
The large y behavior of the scaling function can be 

obtained from the large I behavior of p(g,l) at fixed g. 
Performing an analysis similar to the one leading to (jH 
EJ one finds p(g,l) ~ l~ X ~ 1/tl [h{g)] 2 which, for large g, 
behaves like 

p(g,l)~g-»r 1 - 1 ^ , l,g^oo. (22) 

It follows immediately that for large y, p(y) ~ y~ ~ m as 
announced in Eq. (|6]). This scaling form can be shown to 
be consistent with the large I behavior of Ptime(0 in Eq. 
(gj). We have checked that Eq. is corroborated by 
numerical simulations with a good accuracy for different 
values of \x. 

In this letter, we have investigated the statistical prop- 
erties of the first gap, G n , and the associated time in- 
terval between the two rightmost positions, L„, for a 
one-dimensional random walk in the limit of infinitely 
many time steps. The scaling forms of the limiting 
pdf's p(g,l), Pgap(ff), and ptime(0 for large g and I have 
been obtained, see Eqs. d3j) to (J6j) . Remarkable, unex- 
pected, results are the freezing of the tail of Ptimo(0 to 
Ptime(0 ~ l~ 2 for < /i < 1, and the divergence of the 
average duration (I) — J I Pti mc (l) dl for any Levy index 
< fx < 2. Moreover, while the first moment of the gap 
(g) = J IptimeQ) dg is finite for 1 < /1 < 2, we found that 
it diverges for < fi < 1. Such a divergence is at vari- 
ance with the empirical Bath's law [l6| for earthquakes 
that predicts a finite average gap of magnitude between 
the mainshock and the next largest aftershock. While we 
focused here on the first gap, a recent study of the case 
fx = 2 26| showed that the k th gap, with k large (and 
n — > 00), displays universal fluctuations with a power 
law tail. This naturally raises the question about the 
time between the k th and (k + l) th maximum in the limit 
of large k, which is a challenging problem. 



5 



SNM and GS acknowledge support by ANR grant 
2011-BS04-013-01 WALKMAT and in part by the Indo- 
French Centre for the Promotion of Advanced Research 
under Project 4604-3. 



* majumdar@lptms. u-psud.fr| 

t philippe.mounaix@cpht.polytechnique.fr 

* gregory.schehr@lptms.u-psud.fr 

[1] E. J. Gumbel, Statistics of Extremes, Dover, (1958). 

[2] R. W. Katz, M. P. Parlange and P. Naveau, Adv. Water 

Resour. 25, 1287 (2002). 
[3] P. Embrecht, C. Kluppelberg, T. Mikosh, Modelling 

Extremal Events for Insurance and Finance (Springer), 

Berlin (1997). 

[4] S. N. Majumdar, J.-P. Bouchaud, Quant. Fin. 8, 753 
(2008). 

[5] J.-P. Bouchaud, M. Mezard, J. Phys. A 30, 7997 (1997). 
[6] D. S. Dean, S. N. Majumdar, Phys. Rev. E 64, 046121 
(2001). 

[7] P. Le Doussal and C. Monthus, Physica A 317, 140 
(2003). 

[8] F. Omori, J. Coll. Sci. Imp. Univ. Tokyo 7, 111 (1894). 
[9] T. Utsu, Geophys. Mag. 30, 521 (1961). 
[10] F. billo, R. N. Mantegna, Phys. Rev. E 68, 016119 

(2003) . 

[11] A. M. Petersen, F. Wang, S. Havlin, H. E. Stanley, Phys. 

Rev. E 82, 036114 (2010). 
[12] C. Monthus, P. Le Doussal, Eur. Phys. J. B 41, 535 

(2004) . 

[13] A. G. Pakes, Y. Li, Statist. Probab. Lett. 40, 395 (1998). 
[14] S. Sabhapandit, S. N. Majumdar, Phys. Rev. Lett. 98, 



140201 (2007) 

[15] S. Sabahpandit, S.N. Majumdar, S. Redner, J. Stat. 

Mech. L03001, (2008) 
[16] M. Bath, Techtonophysics 2, 483 (1965). 
[17] D. Vere- Jones, Bull. Seism. Soc. Am. 59, 1535 (1969). 
[18] R. Console, A. M. Lombardi, M. Murru, D. Rhoades, J. 

Geophys. Res. 108, 2128 (2003). 
[19] A. Helmstetter, D. Sornette, Geophys. Res. Lett. 30, 

2069 (2003) and Ref. therein. 
[20] H. A. David, H. N. Nagaraja, Order Statistics (third ed.), 

Wiley, New Jersey (2003). 
[21] W. Yang and Y. Ben-Zion, Geophys. J. Int. 177, 481 

(2009) . 

[22] P. Le Doussal, A. Rosso, K. J. Wiese, Europhys. Lett. 

96, 14005 (2011). 
[23] E. Brunet, B. Derrida, Europhys. Lett. 87, 60010 (2009); 

J. Stat. Phys. 143, 420 (2011). 
[24] N. R. Moloney, K. Ozogany, Z. Racz, Phys. Rev. E 84, 

061101 (2011). 

[25] S. N. Majumdar, R. Ziff, Phys. Rev. Lett. 101, 050601 
(2008). 

[26] G. Schehr, S. N. Majumdar, Phys. Rev. Lett. 108, 040601 
(2012). 

[27] For /i = 1/k, with k £ N, there are some logarithmic 

corrections to the algebraic decay. 
[28] A. Comtet, S. N. Majumdar, J. Stat. Mech. P06013, 

(2005). 

[29] F. Pollaczeck, C. R. Acad. Sci. Paris, 234, 2334 (1952); 

J. Appl. Probab. 12(2), 390 (1975). 
[30] F. Spitzer, Trans. Am. Math. Soc. 82, 323 (1956). 
[31] For a review see S. N. Majumdar, Physica A 389, 4299 

(2010) . 

[32] V. V. Ivanov, Astron. Astrophys. 286, 328 (1994). 
[33] S. N. Majumdar, A. Comtet, R. Ziff, J. Stat. Phys. 122, 
833 (2006). 



